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Abstract 

An  iterative  technicue  and  the  joint  epicenter  determination  tech¬ 
nique  are  utilized  to  compute  improved  locations  for  33  presumed  under¬ 
ground  explosions  in  the  western  Soviet  Union  and  station  corrections  to 
the  Herrin  et  al  (1968)  travel  time  curves  for  98  stations.  Comparison 
of  the  results  of  the  two  techniques  indicates  that  the  iterative  pro¬ 
cedure  produces  improved  location  accuracy  relative  to  the  joint  epi¬ 
center  determination  technique.  Analysis  of  the  station  corrections 
determined  by  this  study  indicate  that  epicenters  located  on  a  worldwide 
basis  with  these  corrections  are  as  reliable,  or  slightly  better  than 
those  determined  using  corrections  from  previous  studies.  For  events  in 
the  western  Soviet  Union,  the  station  corrections  are  superior  to  previ¬ 
ous  estimates.  Comparison  of  the  travel  time  corrections  with  magnitude 
residuals  previously  determined  by  North  (1977)  fails  to  indicate  a 
strong  worldwide  correlation  of  early  arrivals  with  increased  amplitudes 
and  late  arrivals  with  decreased  amplitudes. 
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I .  INTRODUCTION 


A  considerable  effort  has  been  expended  In  the  past  twenty  years  to 
refine  and  explain  the  variations  observed  In  the  travel  time  of  P-waves 
at  teleseismic  distances.  The  early  travel-time  curves  were  based  pri¬ 
marily  on  observations  from  a  large  number  of  earthquakes  at  a  limited 
number  of  seismic  stations.  Later  work  involved  the  use  of  underground 
explosions  for  which  the  locations  and  origin  times  were  precisely  known 
and  recordings  were  made  at  well  distributed  networks  of  stations. 

Utilizing  a  number  of  underground  explosions  and  a  large  suite  of 
well  located  earthquakes,  Herrin  et  al  (1968)  developed  a  revised  travel 
time  curve  which  was  to  represent  an  average  earth.  Many  authors 
(Cleary  and  Hales  (1966);  Engdahl ,  Sindorf  and  Eppley  (1977);  Hales  et 
al  (1968);  Hales  and  Doyle  (1967);  Herrin  and  Taggart  (1968);  Lilwall 
and  Douglas  (1970);  Hales  and  Herrin  (1972)  and  others)  have  examined 
the  computation  of  corrections  to  the  travel  time  curves.  These  correc¬ 
tions  have  been  interpreted  in  various  ways;  relating  them  to  upper  man¬ 
tle  inhomogeneity,  near  station  effects  including  azimuthally  varying 
corrections,  and  also  as  refinements  to  the  actual  average  earth  rela¬ 
tionships.  The  overriding  conclusion  is  that  in  order  to  conduct 
extremely  accurate  epicenter  location  calculations,  specific,  but  vari¬ 
able,  travel  time  corrections  must  be  developed  for  each  station.  When 
a  "worldwide"  average  correction  is  developed  for  a  station,  an  improve¬ 
ment  in  location  capability  is  achieved.  However,  significantly 
improved  location  capabilities  only  come  through  the  development  of  sta¬ 
tion  corrections  for  paths  to  each  geographical  area. 
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It  is  the  purpose  of  this  paper  to  present  the  results  of  an 
analysis  of  travel  time  corrections  to  the  Herrin  et  al  (1968)  curves 
for  a  worldwide  network  of  stations  for  events  in  the  western  Soviet 
Union.  Through  the  use  of  these  corrections,  improved  location  determi¬ 
nation  of  events  in  this  area  will  be  possible.  Two  techniques,  an 
iterative  epicenter  determination  procedure  and  the  Joint  Epicenter 
Determination  (Douglas,  Tilwall  and  Young  (1974))  are  utilized  to 
develop  travel  time  corrections  based  upon  13  presumed  underground 
explosions.  The  two  techniques  are  compared  for  this  application  and 
the  resulting  travel  time  corrections  are  related  to  previous  tabula¬ 
tions  for  other  geographical  regions. 

TI .  TRAVKI,  TIME  CORRECTION  PROCEDURE 

The  procedure  utilized  to  locate  the  presumed  explosions  and  deter¬ 
mine  the  travel  time  corrections  is  quite  similar  to  that  used  by  Herrin 
et  al  (1968)  whereby  residuals  are  developed  from  best  estimates  of 
hypocenters.  Corrections  are  then  developed  from  the  residuals  and 
applied  to  the  data  prior  to  relocating  the  epicenters.  In  this  case, 
measurements  for  epicentral  distances  less  than  20  degrees  and  PKP  dis¬ 
tances  are  rejected.  The  hypocenter  program  (Cannon  (1967))  which  was 
used  in  this  studv,  estimates  the  event  location  through  a  least  squares 
procedure  and  was  set  to  automatically  reject  a  station  when  the  travel 
time  residual  for  the  -'chase  was  greater  than  three  standard  deviations 
from  the  predicted  travel  time  for  the  particular  event.  The  residuals 
for  all  events  were  sorted  bv  station  and  any  phase  more  than  three 


standard  deviations  iron  the  mean  residual  of  the  station  was  also 


rejected.  The  mean  station  residuals  for  the  selected  geographical 
region  were  then  determined  as  the  average  for  each  station.  Correc¬ 
tions  for  these  residuals  were  then  applied  to  the  initial  data  base, 
including  those  phases  which  were  previously  rejected  for  having  large 
residuals.  The  procedure  was  repeated  using  the  corrected  arrival 
times. 

A  data  base  of  31  presumed  underground  explosions  in  the  western 
Soviet  Union  was  selected  from  tabulations  by  Dahlman  and  Israelson 
(1977).  Presumed  explosions  were  chosen  so  that  depth  could  be  con¬ 
strained  to  sea  level  and  thereby  minimize  depth  uncertainty.  The 
events  were  chosen  to  maximize  coverage  of  a  specific  geographical 
region,  namely  the  western  USSR,  thus  minimizing  any  specific  aspect  of 
a  local  region,  and  to  include  events  whose  sizes  would  ensure  that  they 
were  relatively  well  recorded  at  a  large  number  of  stations.  The  dis¬ 
tribution  of  these  events  is  indicated  in  Figure  1.  Available  seismo- 
logical  bulletins  from  the  US  Geological  Survey,  International  Seismo- 
logical  Centre,  VELA  Uniform  Program,  and  several  other  international 
sources  were  searched  for  all  available  P-phase  arrival  times  for  the 
event  set.  This  resulted  in  a  tabulation  of  4,893  reported  phases  at 
835  different  stations.  Only  tabulated  arrival  times  as  reported  by  the 
various  bulletins  were  used  in  this  study.  In  order  to  eliminate  those 
stations  which  detected  a  small  number  of  events,  only  those  which 
detected  at  least  half  of  the  events  were  selected.  This  reduced  the 
number  of  detections  to  2,361  and  the  number  of  stations  to  102.  Four 
additional  stations  were  deleted  when  preliminary  estimations  of  the 


standard  deviation  of  their  travel  *■  ; ,.,u  residuals  were  found  to  he  sig¬ 
nificant1;  larger  than  the  majority  of  the  data  set.  These  stations 
(KJN,  NUR,  SOD,  and  TEH)  were  found  to  also  have  the  smallest  average 
source-to-receiver  distance.  Their  standard  deviations  were,  therefore, 
large  because  of  increased  travel  time  scatter  observed  at  regional  dis¬ 
tances.  As  a  result  of  these  contraints,  the  final  data  set  included 
2,262  detections  at  98  stations.  Less  than  9  percent  of  this  data  set 
was  later  rejected  for  exceeding  the  three  standard  deviation  limits  in 
the  data  reduction  phase  of  the  task.  The  98  station  network  averaged 
63  stations  per  event,  ranging  from  a  minimum  of  25  to  a  maximum  of  87. 

The  iterative  procedure  described  above  was  followed  while  several 
estimates  of  the  convergence  of  the  method  were  being  monitored.  The 
shift  in  the  epicenters  from  one  iteration  to  the  next  was  noted  to  be 
fairly  random  for  the  first  three  iterations  and  then  to  degenerate 
rapidly  to  shifting  all  epicenters  approximately  0.8  kilometers  in  the 
same  direction  during  successive  iterations.  The  mean  of  the  absolute 
values  of  the  travel  time  residual  for  each  iteration  decreased  rapidly 
for  three  iterations  and  then  held  fairly  constant.  Other  criteria  for 
monitoring  the  convergence  indicated  a  similar  relationship  and  it  was 
determined  that  the  third  iteration  results  were  as  accurate  as  the  pro¬ 
cedure  was  capable  of  producing  for  the  given  data  set. 

In  order  to  assess  the  effect  of  a  disproportionate  number  of  sta¬ 
tions  in  North  America  and  Europe,  a  reduced  network  of  33,  averaging  22 
stations  per  event,  was  selected  to  provide  the  best  azimuthal  coverage 
possible  within  the  constraints  of  the  initial  data  base.  Figure  2 
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indicates  the  locations  of  both  the  33  and  98  station  networks.  The 


same  pattern  of  convergence  to  the  third  iteration  was  observed  for  the 
33  station  network  as  with  the  larger  network.  In  most  cases,  the 
resulting  epicenters  were  within  the  95  percent  confidence  ellipses  of 
the  larger  network.  However,  the  areas  of  the  95  percent  confidence 
ellipses  were  approximately  40  percent  larger  in  the  reduced  33  station 
network.  This  variation  in  confidence  ellipses  is  nearly  as  expected 
since  the  square  root  of  the  number  of  stations  is  used  in  computing  the 
ellipses.  Taking  this  effect  into  account,  there  is  no  observable 
difference  in  confidence  region  for  the  two  networks.  This  indicates 
that  the  precision  of  the  procedure  utilized  to  determine  epicenters  and 
estimate  travel  time  residuals  is  fairly  insensitive  to  concentrations 
of  stations  provided  that  the  widest  possible  azimuthal  variation  is 
utilized.  The  results  of  these  two  sets  of  calculations  are  tabulated 
in  Tables  I  and  II. 

It  is  difficult  to  evaluate  the  absolute  accuracy  of  the  epicenters 
which  are  determined  through  this  procedure.  The  confidence  ellipses 
which  are  used  m  estimate  the  location  accuracy  are,  in  actuality, 
measuring  the  fit  of  the  travel  time  data  and  are  not  related  to  the 
actual  locations  of  the  events.  In  order  to  assess  the  true  accuracy  of 
the  procedure,  the  exact  location  of  an  event  must  be  known.  Because 
the  Soviets  do  not  release  specific  details  of  their  underground  nuclear 
testing  program,  we  must  resort  to  circumstantial  evidence.  Of  the  33 
events  utilized  in  this  study,  the  location  of  only  one  is  known  with 
some  degree  of  confidence.  This  event  was  described  by  Marshall  (1972); 
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Table  II 


Travel  Time  Corrections  Resulting  from  Iterative  Procedure 

98  Station  Network  33  Station  Network 


Travel 

Travel 

Time 

Standard 

Time 

Standard 

Correction 

Deviation 

No 

Correction 

Deviation 

No 

Station 

(sec) 

(sec) 

Events 

(sec) 

(sec) 

Events 

ALE 

-0.1 

0.55 

22 

-0.1 

0.50 

22 

ALQ 

-0.5 

0.47 

18 

AVE 

0.0 

0.71 

16 

BER 

-0.2 

0.66 

20 

BLC 

0.4 

0.23 

23 

0.4 

0.22 

23 

BMO 

0.6 

0.53 

32 

0.5 

0.53 

32 

BNG 

1.0 

0.55 

21 

0.7 

0.53 

21 

BNS 

-0.7 

0.53 

16 

BRA 

0.2 

0.82 

12 

BSF 

-0.2 

0.61 

18 

BUH 

-0.6 

0.61 

18 

BUL 

0.7 

0.50 

26 

0.4 

0.52 

26 

CLL 

0.0 

0.32 

30 

CMP 

-2.0 

0.66 

15 

COL 

-0.7 

0.46 

30 

-0.8 

0.44 

30 

COP 

-0.3 

0.50 

17 

CPO 

0.9 

0.42 

23 

0.8 

0.45 

23 

DOU 

-0.5 

0.61 

16 

DUG 

-1.0 

0.47 

21 

DUR 

-0.5 

0.74 

13 

EDM 

-0. 1 

0.26 

26 

-0.1 

0.26 

26 

EKA 

0.1 

0.44 

23 

0.0 

0.49 

23 

ESK 

0.2 

0.54 

16 

EUR 

-0.9 

0.41 

30 

-0.9 

0.40 

30 

FCC 

0.2 

0.51 

18 

FFC 

0.5 

0.51 

25 

FLN 

0.1 

0.49 

26 

0.0 

0.44 

26 

FSJ 

-0.7 

0.36 

16 

FUR 

-1.0 

0.63 

21 

-1.1 

0.56 

20 

GBA 

0.5 

0.51 

21 

0.1 

0.41 

21 

GDH 

0.1 

0.50 

16 

GIL 

-0.7 

0.39 

24 

GMA 

-0.1 

0.42 

16 

GOL 

-0.1 

0.39 

15 

GRF 

-1.1 

0.38 

26 

GRR 

0.3 

0.47 

25 

HFS 

0.9 

0.77 

25 

HYB 

0.9 

0.67 

18 

0.5 

0.54 

18 

IFR 

0.5 

0.68 

20 

0.3 

0.57 

20 

INK 

0.2 

0.43 

21 
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Table  II  (Continued) 


Travel  Time  Corrections  Resulting  from  Iterative  Procedure 

98  Station  Network  33  Station  Network 


Travel 

Travel 

Time 

Standard 

Time 

Standard 

Correct  ion 

Deviation 

No 

Correction 

Deviation 

No 

Station 

(sec) 

(sec) 

Events 

(sec) 

(sec) 

C 

<u 

> 

w 

ISK 

-0.2 

0.48 

13 

JAS 

-0.6 

0.43 

22 

KAS 

-0.6 

0.82 

15 

KBL 

-0.3 

0.64 

18 

-0.5 

0.62 

20 

KBS 

-0.1 

0.47 

17 

-0.2 

0.54 

18 

KDC 

0.6 

0.50 

23 

0.5 

0.50 

23 

KEV 

0.6 

0.64 

27 

0.  7 

0.56 

27 

KHC 

-0.6 

0.  39 

31 

KIC 

1.0 

0.46 

20 

0.7 

0.44 

20 

KIR 

0.9 

0.67 

26 

KJF 

1.1 

0.80 

13 

KON 

0.2 

0.  37 

25 

KRA 

0.0 

0.62 

21 

KRR 

0.6 

0.52 

19 

0.3 

0.53 

19 

KTG 

-1.3 

0.26 

15 

LAO 

-0.8 

0.22 

29 

LBF 

0.4 

0.30 

19 

LJU 

-1.2 

0.41 

19 

LMR 

-0.2 

0.50 

19 

LNS 

0.2 

0.  72 

16 

LON 

-0.2 

0.41 

15 

LOR 

0.4 

0.27 

28 

MAT 

1.1 

0.51 

21 

0.9 

0.40 

21 

MBC 

-0.4 

0.39 

28 

-0.4 

0.34 

28 

MNY 

0.3 

0.55 

16 

MOX 

-0.4 

0.41 

32 

-0.6 

0.35 

32 

NAO 

1.2 

0.  52 

15 

NDI 

0. 1 

0.50 

24 

-0.3 

0.46 

24 

NEW 

-0.2 

0.40 

25 

NIE 

-0.2 

0.80 

16 

NT! 

-0.3 

0.  50 

21 

PMR 

0.6 

0.  39 

26 

PNT 

0.2 

0.36 

21 

PRA 

-0.9 

0.61 

19 

PRU 

-0.7 

0.41 

30 

QUF. 

-0.2 

0.  o3 

16 

-0.7 

0.57 

16 

RES 

0.1 

0.52 

14 

0. 1 

0.50 

14 

RSL 

-0. 1 

0.43 

18 

SCH 

0.6 

0.42 

18 

0.5 

0.40 

18 

SDB 

0.4 

0.  38 

18 

0. 1 

0.32 

18 
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Table  II  (Continued) 


Travel  Time  Corrections  Resulting  from  Iterative  Procedure 

98  Station  Network  33  Station  Network 

Travel  Travel 


Time 

Standard 

Time 

Standard 

Correct  ion 

Deviation 

No 

Correction 

Deviation 

No 

Station 

(sec) 

(sec) 

Events 

(sec) 

(sec) 

Events 

SES 

0.4 

0.41 

27 

SHI 

-0.3 

0.57 

19 

-0.8 

0.54 

19 

SHL 

0.8 

0.63 

15 

0.7 

0.68 

15 

SPF 

-0.3 

0.  35 

19 

SSC 

0. 1 

0.  37 

25 

SSF 

0.3 

0.38 

17 

TCF 

-0.2 

0.  33 

18 

TFO 

-1.2 

0.36 

30 

TRO 

0.6 

0.59 

25 

TUC 

-0.7 

0.63 

15 

-0.8 

0.64 

15 

TUL 

0. 1 

0.46 

25 

UBO 

-0.2 

0.40 

24 

UME 

0.9 

0.57 

22 

UPP 

1.2 

0.66 

26 

VIE 

-1.3 

0.62 

18 

VKA 

-1.2 

0.  57 

20 

VRI 

-0.2 

0.  72 

11 

-0.5 

0.79 

13 

YKC 

0.5 

0.42 

23 

Mean  0.504 

Standard  Deviation  0.137 

Rejected  Stations 


K  JN 

1.33 

13 

NUR 

1  .61 

23 

SOD 

0.99 

31 

TEH 

1.06 

1  1 

1  1 


Dahlman  and  Israelson  ( 1  *) 7 7 ) ;  and  Nordyke  (  1973)  as  being  a  cratering 
experiment  in  a  river  bed  near  49. 9N,  79. OE  on  15  Jan  1965.  As  a  result 
of  the  experiment,  Dahlman  and  Israelson  (1977)  report  that  a  crater  408 
meters  in  diameter  was  formed  and  the  placement  of  the  crater  in  the 
river  bed  caused  the  formation  of  a  reservoir  in  the  river.  A  later 
excavation  of  the  lip  of  the  crater  caused  the  crater  itself  to  also 
fill  with  water  (Nordyke,  1973).  A  dammed  river  and  accompanying  circu¬ 
lar  water-filled  feature  approximately  400-450  meters  in  diameter  can  be 
observed  in  a  I.ANDSAT  photograph  taken  on  15  Aug  1975.  The  geodetic 
coordinates  of  the  crater  are  49.950  degrees  N  by  79.010  degrees  E. 

This  is  approximately  3.1  kilometers  north  of  the  hypocenter  located 
using  98  stations.  If  we  assume  that  this  crater  is  truly  associated 
with  the  15  Jan  63  event,  then  the  apparent  error  is  well  within  the  95 
percent  coverage  ellipse  which  was  estimated  to  he  22.4  by  21.2  kilome¬ 
ters.  Although  this  accuracy  cannot  be  extrapolated  to  the  other  events 
or  even  to  this  event  absolutely  (the  event  conceivably  could  be  associ¬ 
ated  with  the  wrong  feature  on  the  photograph),  it  does  serve  to  lend 
some  confidence  to  the  locations  and  the  procedure  which  was  utilized. 
HI.  HUNT  EPICENTER  DETERMINATION  COMPARISON 

In  order  to  evaluate  the  relative  usefulness  of  the  iterative  hypo- 
enter  location  technique  which  was  described  above  and  the  joint  epi¬ 
center  determination  procedure,  Douglas,  Lilwall  and  Young  (1974),  a 
reduced  data  base  was  constructed  from  that  used  in  the  iterative  pro¬ 
cedure.  Since  the  joint  epicenter  determination  procedure  computer  pro¬ 
gram,  which  was  obtained  from  the  Seismic  Data  Analysis  Center,  did  not 

12 

l 


have  provisions  for  eliminating  erroneous  data  points,  only  those  P-wave 


phases  which  were  utilized  in  the  final  iteration  of  the  Iterative  tech¬ 
nique  were  included.  Both  the  98  and  33  station  networks  were 
evaluated.  The  results  of  this  analysis  are  tabulated  in  Tables  III  and 
IV. 

A  comparison  of  the  results  of  the  two  techniques  indicates  a  slight 
reduction  in  the  estimated  95  percent  coverage  limits  of  the  joint  epi¬ 
center  determination  procedure.  This  is  to  be  expected  as  pointed  out 
previously  by  Ahner,  Blandford  and  Shumway  (1971),  since  the  joint  epi¬ 
center  procedure  is  based  upon  minimizing  the  error  of  the  travel  times 
of  the  entire  data  set  through  least  squares  procedures  whereas  the 
iterative  reduces  the  error  for  each  event  individually.  It  should  also 
be  noted  that  the  98  station  joint  epicenter  determination  epicenters 
are  shifted  an  average  of  11.6  kilometers  to  the  northwest  of  the  itera¬ 
tive  procedure  epicenters.  The  variation  in  this  shift  is  remarkahlv 
low  as  If  most  epicenters  were  shifted  as  a  group.  A  similar  relation¬ 
ship  exists  with  the  33  station  network  epicenters  as  determined  through 
the  Joint  epicenter  determination  procedure.  The  95  percent  confiderce 
limits  are  slightly  larger  for  the  33  station  network  but  the  shifts 
from  the  iteratively  determined  epicenters  averages  only  7.7  kilometers 
to  the  northwest  with  a  slightly  larger  scatter.  The  98  station  joint 
epicenter  determination  site  for  the  15  .Jan  65  event  is  approximacpi v 
12.5  kilometers  northwest  of  the  LANDSAT  crater.  The  33  station  loint 
epicenter  determination  site  is  approximately  12.2  kilometers  west  of 
the  LANDSAT  crater.  Figure  3  presents  a  map  of  the  locations  determined 
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TABLE  IV 


Travel  Time  Corrections  Resulting  From  Joint  Epicenter  Determination 
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98  STATION  NETWORK  33  STATION  NETWORK 


Travel 

Travel 

Time 

Standard 

Time 

Standard 

Station 

Correction 

Deviation 

No 

Correction 

Deviation 

No 

(sec) 

(sec) 

Events 

(sec) 

(sec) 

Events 

ALE 

0.2 

0.25 

22 

-0.1 

0.33 

22 

ALQ 

-0.3 

0.25 

18 

AVE 

-0.2 

0.26 

16 

BER 

-0.6 

0.25 

20 

BLC 

0.7 

0.23 

23 

0.4 

0.30 

23 

BMO 

0.9 

0.20 

32 

0.5 

0.24 

32 

BNG 

1.1 

0.27 

21 

0.6 

0.30 

21 

8NS 

-1.0 

0.27 

16 

BRA 

-0.3 

0.31 

12 

BSF 

-0.6 

0.25 

18 

BUH 

-1.1 

0.25 

18 

BUL 

1.0 

0.28 

26 

0.3 

0.28 

26 

CLL 

-0.4 

0.21 

30 

CMP 

-2.4 

0.29 

15 

COL 

-0.1 

0.24 

30 

-0.6 

0.28 

30 

COP 

-0.6 

0.26 

17 

CPC 

1.2 

0.22 

23 

0.8 

0.25 

23 

DOU 

-0.8 

0.27 

16 

DUG 

-0.8 

0.24 

21 

mm 

-0.9 

0.29 

13 

EDM 

0.3 

0.22 

26 

-0.1 

0.26 

26 

EKA 

-0.2 

0.23 

23 

-0.3 

0.34 

23 

ESK 

-0.2 

0.26 

16 

EUR 

-0.6 

0.21 

30 

-1.1 

0.24 

30 

FCC 

0.5 

0.25 

18 

FFC 

0.9 

0.22 

25 

FLN 

-0.2 

0.22 

26 

-0.3 

0.32 

26 

FSJ 

-0.3 

0.28 

16 

FUR 

-1.4 

0.24 

21 

-1.5 

0.34 

21 

GBA 

1.6 

0.40 

21 

0.6 

0.44 

21 

GDH 

0.2 

0.27 

16 

GIL 

-0.1 

0.25 

24 

GMA 

0.5 

0.30 

16 

GOL 

0.0 

0.27 

15 

GRF 

-1.5 

0.22 

26 

GRR 

-0.1 

0.22 

25 

KFS 

0.6 

0.24 

25 

KYB 

2.1 

0.42 

18 

1.0 

0.46 

18 

IFF. 

0.0 

0,24 

20 

-0.2 

0.30 

20 

15 


TABLE  IV  (Continued) 


Trivel  Time  Corrections  Resulting  From  Joint  Epicenter  Determination 


98  STATION  NETWORK  33  STATION  NETWORK 


Travel 

Travel 

T  ime 

Standard 

Time 

Standard 

Correct  ion 

Deviation 

No 

Correction 

Deviation 

No 

Station 

(sec) 

(sec) 

Events 

(sec) 

(sec) 

Event 

INK 

0.7 

0.26 

21 

ISK 

-0.6 

0.32 

13 

JAS 

-0. 1 

0.24 

11 

KAS 

-0.9 

0.  32 

15 

KBL 

0.6 

0.44 

18 

-0.5 

0.51 

18 

KBS 

0.0 

0.  28 

17 

-0.2 

0.37 

17 

KDC 

1.3 

0.26 

23 

0.8 

0.29 

23 

KEV 

0.5 

0.24 

27 

0.4 

0.37 

27 

KHC 

-1.1 

0.21 

31 

KIC 

1.0 

0.25 

20 

0.6 

0.28 

20 

KIR 

0.8 

0.24 

26 

KJF 

0.9 

0.31 

13 

KON 

-0.2 

0.23 

25 

KRA 

-0.4 

0.25 

21 

KRR 

0.9 

0.30 

19 

0.2 

0.30 

19 

KTG 

-1.3 

0.28 

15 

LAO 

-0.4 

0.21 

29 

LBF 

-0.1 

0.25 

19 

LJU 

-1.6 

0.25 

19 

LMR 

-0.5 

0.25 

19 

LNS 

-0.3 

0.27 

16 

LON 

0.2 

0.28 

15 

LOR 

0.0 

0.21 

28 

MAT 

2.5 

0.36 

21 

1.6 

0.40 

21 

M3C 

0.1 

0.23 

28 

-0.3 

0.30 

28 

MNY 

0.0 

0.27 

16 

MOX 

-0.9 

0.20 

32 

-1.0 

0.32 

32 

NAO 

0.8 

0.28 

15 

ndi 

1.3 

0.42 

24 

0.2 

0.48 

24 

NEW 

0.3 

0.23 

25 

NIE 

-0.7 

0.28 

16 

NT  I 

-0. 1 

0.24 

21 

PMR 

1.3 

0.25 

26 

PNT 

0.6 

0.24 

21 

PRA 

-1.3 

0.25 

19 

PRU 

-1.1 

0.21 

30 

Q'JE 

0.5 

0.44 

16 

-0.6 

0.50 

16 

RES 

0.4 

0.29 

14 

0.1 

0.35 

14 

RSL 

-0.6 

0.25 

18 

16 


TABLE  IV  (Continued) 

Travel  Time  Corrections  Resulting  From  Joint  Epicenter  Determination 


Station 


SCH 

SDB 

SES 

SHI 

SHL 

SPF 

SSC 

SSF 

TCF 

TFO 

TRO 

TUC 

TUL 

UBO 

UHE 

UPP 

VTE 

VKA 

VRI 

YKC 


98  STATION  NETWORK  33  STATION  NETWORK 


Travel 

Time 

Standard 

Correction 

Deviation 

(sec) 

(sec) 

0.6 

0.25 

0.5 

0.29 

0.8 

0.22 

0.0 

0.38 

2.3 

0.45 

-0.6 

0.25 

-0.3 

0.22 

0.0 

0.21 

-0.6 

0.25 

-0.9 

0.20 

0.4 

0.24 

-0.3 

0.27 

0.5 

0.21 

0.0 

0.22 

0.6 

0.25 

0.8 

0.23 

-1.7 

0.26 

-1.7 

0.25 

-0.8 

0.33 

0.9 

0.24 

Travel 

Time 

No 

Correction 

Events 

(sec) 

18 

0.4 

18 

-0.1 

27 

19 

-0.8 

15 

1.1 

19 

25 

27 

18 

30 

25 

15 

-0.8 

25 

24 

22 

26 

18 

20 

11 

-1.0 

23 

Standard 

Deviation 

No 

(sec) 

Events 

0.31 

18 

0.29 

18 

0.42 

19 

0.49 

15 

0.29  15 


0.42  11 


17 


I 


N  Latitude 


for  this  event.  In  both  of  these  eases,  the  location  error  exceeds  the 


95  percent  confidence  limits  hy  a  factor  of  two,  if  we  assume  that  the 
crater  and  event  are  truly  associated.  This  would  indicate  that  the 
iterative  procedures  have  produced  more  accurate  locations  with  errors 
within  the  stated  bounds. 

IV.  COMPARISON  TO  PREVIOUS  STUDIES 

The  comparison  of  travel  time  corrections  derived  through  various 
procedures  for  different  travel  time  curves  and  geographical  regions  is 
of  interest  for  several  reasons.  It  is  primarily  of  importance  to  esti¬ 
mate  the  variations  which  may  be  encountered  between  different  geograph¬ 
ical  areas  and  in  analyzing  different  travel  time  curves.  Through  the 
analysis  of  these  types  of  variations,  it  may  be  possible  to  estimate 
the  size  of  an  area  over  which  the  travel  time  corrections  remain  valid. 

Comparisons  of  the  residuals  determined  through  the  iterative  pro¬ 
cedure  for  the  98  station  network  have  been  made  with  summaries  by 
Cleary  and  Hales  (1966),  Lilwall  and  Douglas  (1970),  Sengupta  and  Julian 
(1976),  and  Masso,  Savino,  and  Bache  (1978).  As  noted  in  Figures  4  and 
5,  the  Lilwall  and  Douglas  (1970)  and  Sengupta  and  Julian  (1976)  residu¬ 
als  appear  to  have  the  largest  scatter  of  those  examined.  The  Lilwall 
and  Douglas  (1970)  study,  using  the  joint  epicenter  determination 
method,  also  included  azimuthal  variations  in  the  corrections  along  with 
corrections  to  the  travel  time  relationship.  It  is,  therefore,  not  com¬ 
pletely  relevant  to  consider  their  corrections  the  poorest  based  upon 
the  scatter  on  the  plot. 
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It 


Lilwall  8  Douglas  Residual  (sec) 

(early)  (late) 


(early)  (late) 

98-sta  Iterative  Residual  (sec) 


figure  4.  Comparison  of  Travel  Time 
Residuals  Between  Lilwall 
and  Douglas  and  This  Study 


The  scatter  of  the  residuals  observed  in  this  study  with  those  by 
Cleary  and  Hales  (1966),  Figure  6,  are  more  comparable.  Sengupta  and 
Julian  (1976)  rates  their  residuals,  which  were  based  upon  a  worldwide 
distribution  of  deep  earthquakes  analyzed  through  an  iterative  procedure 
very  similar  to  the  one  used  here,  significantly  better  than  other  stu¬ 
dies  since  the  mean  of  the  standard  deviations  of  their  residuals  was 
0.6  seconds,  much  smaller  than  other  studies.  The  mean  of  the  standard 
deviations  of  the  residuals  in  this  study  is  0.5  seconds  and  by  this 
criteria  is,  therefore,  an  improvement  over  the  previous  studies,  espe¬ 
cially  when  applied  to  the  western  Soviet  Union  source  region. 

Another  recent  study  by  Masso,  Savino  and  Bache  (1978)  analyzed  the 
International  Seismological  Centre  bulletins  from  1964  through  1970  for 
events  with  mb  >  5.0  as  observed  at  524  worldwide  stations.  The 
scatter,  as  observed  in  Figure  7,  between  their  set  of  residuals  for  a 
worldwide  event  set  is  slightly  smaller  in  comparison  to  the  results  of 
this  study  than  any  other  examined.  It  is,  therefore,  conceivable  that 
the  corrections  determined  in  this  study  may  produce  quite  accurate  epi- 
central  estimates  for  events  outside  the  geographical  spread  of  events 
used  in  the  study. 

A  concept  which  has  been  discussed  in  some  detail  in  recent  years 
relates  to  the  relationship  of  travel  time  and  magnitude  residuals.  It 
is  commonly  held  that  early  arriving  signals  correlate  with  high  ampli¬ 
tudes  and  late  arriving  signals  correlate  with  more  attenuated  signals. 
The  classic  example  of  this  is  North  America  where  signals  in  the 
southwestern  United  States  are  late  and  small  relative  to  signals 
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Masso,  Savino  8  Bache  Residual  (sec) 
(early)  (late) 


recorded  In  the  east.  Similar  correlations  to  high  heat  flow  with  the 


higher  attenuation  and  variations  in  Pn  velocities  appear  very  convinc¬ 
ing.  Analysis  of  the  travel  time  residuals  determined  hy  this  study 
with  the  magnitude  residuals  determined  by  North  (1977),  corrected  to  a 
granite  crust,  as  outlined  by  Der  et  al  (1978),  does  not  result  in  a 
similar  conclusion,  as  seen  in  Figure  8.  It  may  be  that  different  geo¬ 
graphical  regions  behave  differently  with  respect  to  these  factors.  The 
commonly  held  relationships  are  observed  in  North  America  where  they 
were  first  put  forth.  However,  European  stations  are  observed  to  be  all 
on  the  high  side  of  the  magnitude  distribution  with  travel  time  residu¬ 
als  fairly  evenly  balanced  between  early  and  late.  Africa  is  distinc¬ 
tively  early  and  low,  contrary  to  the  established  relationship.  This 
anomalous  region  in  Africa  was  previously  pointed  out  by  Masso,  Savino 
and  Raehe  (1978).  This  study  indicates  that  correlations,  on  a  world¬ 
wide  basis,  between  travel  time  and  magnitude  residuals  are  not  as  con¬ 
clusive  as  magnitude  residual  versus  Pn  velocity  and  may  actually 
represent  two  different,  but  related,  effects. 

V.  CONCLUSIONS 

A  set  of  travel  time  correctons  have  been  developed  using  33  events 
in  the  western  Soviet  Union.  The  use  of  these  corrections  at  the  98 
stations  examined  should  make  possible  improved  locations  of  events  in 
this  region.  The  improved  locations  of  the  33  events  in  this  study  are 
also  included. 

In  developing  the  travel  time  corrections,  an  iterative  hypocenter 
location  procedure  and  the  joint  epicenter  determination  procedure  were 
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Figure  8.  Travel  Time  Residuals 

versus  Magnitude  Residuals 


evaluated.  From  the  analysis  of  one  known  event  location.  It  was  deter¬ 
mined  that  the  Iterative  technique  provided  more  accurate  event  loca¬ 
tions  and  travel  time  corrections.  It  was  also  noted  that  the  reduction 
of  the  98  station  network  to  13  stations  to  provide  a  more  uniform  sta¬ 
tion  pattern  did  not  improve  the  results  of  the  procedures.  It  is, 
therefore,  to  be  concluded  that  If  an  adequate  azimuthal  variation  is 
included  in  the  network,  concentrations  of  stations  in  certain  areas  do 
not  significantly  bias  the  results  of  the  iterative  procedures. 

An  analysis  of  the  travel  time  corrections  derived  in  this  study 
with  those  of  previous  studies  indicates  that  these  results  are  as  good, 
or  better,  than  previously  reported  worldwide  travel  time  corrections 
and  that  for  the  geographical  region  investigated,  they  are  superior. 
Comparison  of  the  travel  time  residuals  with  previously  determined  mag¬ 
nitude  residuals  fails  to  indicate  a  strong  worldwide  correlation  of  the 
relationship  of  early  arriving  signals  with  increased  amplitudes  and 
late  arriving  signals  with  reduced  amplitudes  as  previously  observed  in 
North  America. 
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